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Thickness dependence of magnetic properties of (Ga.Mn)As 
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We report on a monotonic reduction of Curie temperature in dilute ferromagnetic semiconductor (Ga,Mn)As 
upon a well controlled chemical-etching/oxidizing thinning from 15 nm down to complete removal of the ferro- 
magnetic response. The effect already starts at the very beginning of the thinning process and is accompanied 
by the spin reorientation transition of the in-plane uniaxial anisotropy. We postulate that a negative gradient 
along the growth direction of self-compensating defects (Mn interstitial) and the presence of surface donor 
traps gives quantitative account on these effects within the p-d mean field Zener model with adequate mod- 
ifications to take a nonuniform distribution of holes and Mn cations into account. The described here effects 
are of practical importance for employing thin and ultrathin layers of (Ga,Mn)As or relative compounds in 
concept spintronics devices, like resonant tunneling devices in particular. 
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Dilute ferromagnetic semiconductors (DFS), such as 
(Ga,Mn)As, are extensively studied in search for new 
spintronic phenomena and towards potential applications 
in memory and information processing technologies, 1 
Among them, the most application-promising is the 
isothermal control of the magnetic phased and magnetic 
anisotropy change 4 by the externally applied electric 
field, which recently led to application-viable demonstra- 
tion of the electric field induced magnetization switching 
in sub-nanometer thin Fe-Co layers^ Thus, there is a gen- 
eral interest in studying various (new) device configura- 
tions in which ultrathin, less than few nanometers, semi- 
conducting (III,Mn)V layers are incorporated. Therefore 
it is timely to provide experimental information how mi- 
cromagnetism in such thin DFS compares with that of 
thicker layers for which a great deal of information has 
already been acquired^— 

In this letter we investigate changes of Curie temper- 
ature (Tq) and magnetic anisotropy associated with sys- 
tematic thinning of d = 15 nm thick as-grown (Ga,Mn)As 
layers till the complete loss of the ferromagnetic signa- 
tures. The magnitude of Tq is as high as ~ 110 K, which 
points to high sample quality, in particular, to a low den- 
sity of antisite compensating donors^ We establish the 
presence of two thickness regimes, for which the evolu- 
tion of magnetic properties is determined by two distinct 
mechanisms. The first is associated with the presence of 
a gradient in the concentration of interstitial Mn ions. 
The second, operating at d < 4 nm, is brought about 
by surface defects pining the Fermi energy in the mid 
gap region of GaAs, and thus depleting holes^^ Our 
results demonstrate, therefore, why low temperature an- 
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nealing not only reduces the concentration of intersti- 
tial Mn but also "homogenizes" magnetic properties of 
(Ga,Mn)As, as observed in neutron studies* 1 - 1 - Further- 
more, the experimentally determined dependence Tc(d) 
in the second regime substantiates the previous theoret- 
ical model of ferromagnetism in interfacial space charge 
layers of (Ga,Mn)Asi2ii2 We also evidence a thickness in- 
duced 90° rotation ([110] -O- [110]) of an in-plane uniaxial 
easy axis [a spin reorientation transition (SRT)]. 

Two d = 15 nm (Ga,Mn)As films (A and B) have been 
deposited at 200° C by low-temperature (LT) molecu- 
lar beam epitaxy on (GaAs) on (100) substrates buffered 
by 500 nm thick LT-GaAs with a use of arsenic valved 
cracker effusion cell. The concentration of substitutional 
Mn at Ga sites x su b — 6% is estimated by the growth 
rate increase of (Ga,Mn)As in comparison to the GaAs 
buffer.— The samples have not undergone any post- 
growth heat treatment. The thickness dependent data 
are obtained by the controlled thinning of the magnetic 
films via sequential open air oxidation of the superficial 
part of the layer which had its native oxide removed 
beforehand by 30 s dipping in concentrated (as 30%) 
HCL 14 ' 15 As this native oxide restores on an expense of 
the top-most part of the semiconductor film, multiple 
repetition of this method allows for very fine and uni- 
form thinning of even macroscopically large areas. The 
thickness of the reformed oxide depends on the oxida- 
tion time^ and it takes 26 etching-oxidation steps to 
completely remove ferromagnetic signal from the layer A 
which oxidized on open air for 6-9 hours between etching 
and magnetic measurements but only 13 steps are needed 
for the control sample B which oxidized about 24 hours. 
As the observed reduction of Tq to thinning is essentially 
the same for both samples, we narrow this report to the 
main sample A only. 

Atomic force microscopy (AFM) images reassure us 
that the whole process proceeds in a truly planar fash- 
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FIG. 1. (Color online) AFM images for two pieces from the 
same wafer (a) AFM image for as grown sample (b) AFM 
image for 26 times etched sample. 



ion. As indicated in Fig. [T]even after 26 etchings we do 
not find any macroscopic or submicrometer sized features 
suggestive a presence of a statistically relevant number of 
(Ga,Mn)As islands left on the GaAs substrate. Neither 
do we observe any unidirectional surface undulation re- 
cently reported for similar (Ga,Mn)As layers*^ Finally, 
we note that it is indeed oxygen needed to conduct the 
process, as no changes of signal are recorded in time do- 
main during lengthy magnetic measurements which are 
done in a chemically inert helium atmosphere. 

Magnetometry measurements are carried out on a 
home made superconducting quantum interference device 
(SQUID) magnetometer operating down to 5 K and up 
to 5 kOe. Routinely field cooled (FC) measurements at 
the magnetic field H = 1 kOe are performed to assess 
the temperature dependence of the spontaneous moment 
m(T) and to establish the magnitude of the saturation 
moment mg = m(l kOe, 5 K). We find such an assess- 
ment of ms quite satisfactory, as prior to the whole exper- 
iment we checked at 5 K that the moment of the layers 
saturates already below 1 kOe for both major in-plane 
crystallographic directions, [110] and [110]. To establish 
Tq of the layer every FC measurement is followed by 
a thermoremnant (TRM) one carried out on increasing 
temperature at H = until the remnant moment van- 
ishes completely. This set of measurements is performed 
twice for these two in-plane directions, and is repeated 
after every step of thinning. We follow the experimental 
code described recently in Ref. [I?]- 

Figure (2^ exemplifies m(T) collected at every 5th step 
of thinning which evidence a gradual weakening of the 
ferromagnetic response in terms of reduction of both the 
magnitude of magnetic moment (and so of its saturation 
value) and its onset temperature (~ Tq). As in the car- 
rier mediated ferromagnetism^ ms depends primarily on 
the number of these Mn cations which are bound ferro- 
magnetically by holes it is straightforward to assume that 
the drop of ms is caused essentially by the reduction of 
the volume of the layer. Therefore, and on the account of 
the AFM studies, we solely assign this drop of m to the 
step-by-step reduction of the (Ga,Mn)As thickness and 
employ ms to size the thickness of the layer after each 
thinning stage. 

In Fig. ^jp we plot the values of ms versus the num- 
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FIG. 2. (Color online) (a) Temperature dependence of the 
field cooled moment m-pc of layer A measured after every 5th 
stage of the thinning process, (b) Reduction of the saturation 
moment [taken as rriFc(5 K)] versus number of the thinning 
step. Dashed line defines the thinning efficiency 77. Calculated 
layer thickness of the remaining material, d n = 15 — my, is 
given on the right j/-axis. 



ber of the thinning step n. The observed there remark- 
ably linear dependence allows us to establish the aver- 
age thickness loss per each thinning step (77 ~ 0.55 nm), 
which is the thickness of the native oxide forming on the 
free surface of the (Ga,Mn)As between etching and begin- 
ning of measurement, and calculate (Ga,Mn)As thickness 
d n = 15 — nr] nm left after each stage of the thinning pro- 
cess. We note in parenthesis that the same experimental 
procedure gives larger 77 ~ 1 nm for sample B, and, given 
the different time of oxidation, both values agree with 
those reported in Refs. [HI and[l9l. 

Bullets in Fig. [3^ illustrate the main effect observed in 
this study: a continuous decrease of Tq in response to the 
reduction of the layer thickness, and, quite remarkably, 
that the whole process is already effective at the very first 
steps of thinning, that is for d ~ 15 nm. We further note, 
that the whole observed 107 — > 43 K drop of Tq consists 
of two parts suggesting that two different mechanisms 
conspire. According to the previous findings we ex- 
pect that the fast drop of Tq for the lowest thicknesses 
results from a significant reduction of hole concentration 
p due to the depletion zones present at both boundaries 
of the layer. However, in order to account for the initial 
reduction of Tq (d = 15 — >• 5 nm) we need to assume 
an existence of a positive gradient (counting along z) of 
either x c s or p. 

Within the framework of the p-d Zener model, Tq 
depends independently on x e g and p, so an adequately 
strong positive gradient of any of these two quantities 
would explain this behavior. However, we actually need 
a simultaneous reduction of both x c g and p to take 
place on thinning in order to explain simultaneously oc- 
curring SRT, which is heralded by [110] — > [TlO] rota- 
tion of the magnetic easy axis taking place at constant 
temperature^ST— as exemplified for T = 94 K in panels 
bl-b3 of Fig. [3] We underline here that the presence of 
SRT calls for the determination of Tq initially from [110] 
projection of TRM (for d > 10 nm) and from the [110] 
projection afterwards, as it is differentiated in Fig. [3^ 
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FIG. 3. (Color online) (a/top) (Full points) experimentally 
determined values and (dashed and solid lines) results of mod- 
eling of Curie temperature Tc in 15 nm (Ga,Mn)As as a func- 
tion of decreasing layer thickness d. (a/bottom) Hole profiles 
for various values of d calculated by solving the Poisson equa- 
tion with the negative gradient of the concentration of Mn 
at interstitial sites and electrostatic effects present at the top 
and bottom boundaries of the layer (see the text). (bl-b3) 
Squares - [110], bullets - [TlO] projection of the thermorem- 
nant moment at elevated temperatures recorded for: (bl) as 
grown layer, (b2) after 5th step, and (b3) after 10th step of 
thinning. The corresponding layer thickness is quoted. The 
vertical dashed line marks an arbitrary selected temperature 
T = 94 K at which thickness-driven spin reorientation transi- 
tion is observed. Arrows indicate Tc for these cases. 

by use of two different colors and symbols to mark the 
experimental values of Tq . 

We now show that the presence of a negative gradi- 
ent of x\ in the layer is a sufficient requirement to ob- 
tain the requested changes of x e g and p on z and we use 
the p-d Zener model of ferromagnetism in DFS to give 
a numerical account on the scale of the effect. Firstly, 
we fix x su b as depth independent according to the sta- 
bility of the reflection high-energy electron diffraction 
pattern recorded during the growth. We assume then, 
for the sake of simplicity, a linear decrease x\ on z, 
x\{z) = £t(15) + (15 — z)b, b > 0, and calculate accord- 
ingly x cS (z) and p(z) using Fig. 8 of Ref. 7 to evaluate 
xi(15). Now with the problem reduced to just only two 
free parameters, namely: x su h and the gradient magni- 
tude 6, their values are obtained by fitting calculated 
Tc(d) to the experimental data using the following for- 
mulae developed to describe Tq of thin and nonuniform 
layers of (Ga,Mn)As&±2 

T c = J dzT c D \p(z),x eS (z)\ J dz P ^, (1) 

where, due to a rather short phase coherence of holes at 
these temperatures L<f, « 1 nm^ z runs from to d 
in the intervals of 1 nm and the maximal obtained value 
is assign to current Tc(d)££ Here T^P is the Curie tem- 
perature calculated within the conventional p-d Zener 
model 18 at given z, p s = J dzp(z) is the sheet hole den- 



sity. Indeed, as indicated in Fig. |3ji/top by the dashed 
line, the presented above simple model reproduces the 
data remarkably well. The fit yields x su b = 5.7%, the 
value consistent with that established from the layer 
growth rate, and b — 0.023 %/nm, or xj drops from 2.15 
at the beginning of the growth to 1.8% at the end. Impor- 
tantly, the total Mn concentration changes along these 
15 nm rather marginally from 7.8 to 7.45% (a relative 
drop by merely ~5%), what is currently beyond the res- 
olution of even the most advanced direct atomic concen- 
tration profilers like secondary ion mass spectroscopy and 
three-dimensional atomic probed Conversely, we can 
state that, despite being tedious, this method provides 
the most accurate (indirect) assessment of the depth de- 
pendence of the two most numerous Mn species in very 
thin (Ga,Mn)As layers. 

At the final stage we add electrostatic effects due to 
the presence of antisite Asc a donors of concentration Nd 
in the LT-GaAs buffer adjacent to (Ga,Mn)As channel 
and donor-like traps at the free (Ga,Mn)As surface. We 
model the traps by introducing an ever-present topmost 
1 nm region of the layer containing Nj donors. We pin 
the Fermi energy at the midgap of GaAs substrate re- 
siding 50 nm below the bottom of our layer. Then we 
solve the Poisson equation within nextnano 3 package 2 ^ 
for the established already distribution of interstitial dou- 
ble donors x\(z) and look for a set of Nj and TVd that 
reproduces the experimental Tc{d). 

Figured shows the hole distribution profiles for var- 
ious thicknesses (bottom part) and (top part) the cor- 
responding Tq values calculated at particular channel 
thickness (thick line) for Nj = 2.8 x 10 20 cm~ 3 and 
Ns = 5 x 10 19 cm -3 . Although this is not a rigorous 
fit and the results depend to some extent on the mag- 
nitudes of these two adjustable concentrations, we can 
conclude that for generally similar values of Nj and JVd 
to those reported previously by some of us 3 -^ the pre- 
sented here model describe quantitatively the magnitude 
of Tc changes on thinning. 

We are now in a position to address the question why 
there exists the Mni gradient in (Ga,Mn)As, and perhaps 
other (III,Mn)V, at the first place. We argue that this 
is this dense and narrow pocket of the surface donor- 
like states that influences how Mni, the by far domi- 
nating self-compensating defect in (Ga,Mn)As, get dis- 
tributed during the growth. We recall here that both 
the traps and Mni are having a similar densities and are 
positively charged, so they repel each other. But as the 
traps are fixed to the surface and Mni are quite mobile at 
the growth temperatures (which is why the LT annealing 
does work in these compounds), so the Mni are being con- 
stantly pushed back from growth front towards the sub- 
strate. It is beyond the scope of this letter to provide with 
a selfconsistent solution of the full thermodynamics of the 
process, but we do see that the combined push backwards 
by the surface positive charge and simultaneous push for- 
ward from the already accumulated Mni at the deeper 
parts of the already grown layer may/should result in an 
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FIG. 4. (Color online) Literature data on Tc (given in units 
of the effective Mn contraction, x e g) of (Ga,Mn) As plotted as 
a function of thickness d. Open symbols - as grown, bullets - 
after annealing, diamonds - present data. Where not explic- 
itly stated, the x e g values have been calculated according to 
Ref. 0- Essentially the same picture is obtained if Tc alone 
is plotted. 2 ® The thick dashed line is a guide for the eye and 
indicates the main trend of improvement of Tc on lowering d. 



equilibrium negative gradient of Mni, which gets imme- 
diately frozen down once the growth is completed, tem- 
perature reduced and Mni immobilized. Therefore, we 
believe that apart from many possible technical-related 
reasons the proposed here mechanism should be effec- 
tive in all (Ga,Mn)As layers with the gradient coefficient 
b being a decreasing function of the thickness, except, 
perhaps of hydrogen-codoped layers where virtually no 
Mni are expected to form. Similarly, in materials like 
(Ga,Mn)Sb where the surface states are filled by the band 
states there should be no 'push-back' effect and so both 
x\ and p should stay constant along the depth of the 
layer, unless a variation of growth parameters tells the 
system otherwise. 

It is interesting to compare our results to previously 
determined values of Tq in films of various thicknesses d. 
According to data collected in Fig. 2J annealed samples 
tend to show increasing Tc values when d decreases down 
to 20 nm. This trend reflects a decreasing efficiency of 
the low temperature annealing with the layer thickness, 
the effect expected within the model of interstitial Mn 
diffusion^ As already discussed, our results reveals an 
opposite trend in as-grown samples. 

In conclusion, we have performed systematic studies 
how Tc in (Ga,Mn)As depends on thickness, finding that 
the observed nearly 60% drop in the magnitude of Tc 
begins already in 15 nm layers where no depletion due 
to surface donor defects is expected. We have assigned 
this effect to a build-in negative gradient of the Mn in- 
terstitials' concentration and successfully reproduced the 
experimental finding in the frame of the adequately mod- 
ified p-d Zener model of ferromagnetism to the case of 
nonuniform hole and Mn distributions. This result is ex- 
pected to help to understand the properties and behavior 



of various spintronic devices which rely they functional- 
ities on thin or even ultra thin layers of (Ga,Mn)As and 
possibly other DFS. 
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FIG. 5. Literature data on Tc of (Ga,Mn)As plotted as a 
function of thickness d. Red numbers - as grown, bold brown 
numbers - after annealing, diamonds - present data. The thick 
dashed line is a guide for the eye and indicates the main trend 
of improvement of Tc on lowering d. 

Figure [5] shows the same Tc values as those in the 
Fig. 4 of the main paper, but plotted without normal- 
ization with respect to the effective Mn concentration 
x e s = %sub — x\{z), where x su b and x\ are concentrations 
of Mn ions on Ga and interstitial sites, respectively. 
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